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INTRODUCTION 


It has been demonstrated by Wdowczyk and Wolfendale 
(1984) that the energy spectrum of the secondary pions pro- 
duced in p - p interactions with energies up to a few times 
10™ eV can be described well by a scaling violation formula 
of the form 


P(E n )=(E o /205) q((E n /E o ) {eJz05) ) 


oC, k(E o , 205)| 


( 1 ) 


E, 


where the parameter o( , which is a measure of the scaling vio- 
lation, is gradually increasing with energy. The coefficient 
k(E ,205) : =K fT (E )/k( 205) describes the variation of inelasti- 
city coefficient for charged pion production. 

The energy spectra of the others secondary particles 
can be obtained by similar transformations. 

For purpose of numerical calculations the function q(x) 
is usually taken from Fermi lab data at 205 GeV and approxi- 
mated by the following formula (for tt”) 


q(x) = . 692 ( 1-x)^*^ exp(2.185x) + ,807(l-x)^* 2 exp( 2 . 968 x)| 


,( 2 ) 


which is obtained by integration over p of the more general 
expression fitting simultaneously the x t and p, variation for 
the same Fermi lab data at 205 GeV ( Kafka etVi 1977 ). The 
expression has the form 


q(x,P t )=9.624( 1-x)^* 2 exp( - 5.963 exp(-2.968x) ) 

+8.744( l-x)^*^ exp(-6.318 p^ exp(-2, 185x) ) 


( 3 ) 


As a reasonable choice for </. and k variation it is taken 


o( =.815 ((E o / 205) 0,135 -1) ln(E Q /205), 
o(=.072 + .0163 ln(E o /205), 


for E q ^ 7.7 10 GeV 
for E q > 7.7 10 4 GeV. (4) 
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k Trtot^ E o» 205 ) !!! ^ E o / 205) “* 042 » k TT±( E o» 205 ) =: ( E O / 2 * * 0 5 )’ #06 * (5) 

It is also assumed that formula (l) and (2) can be used for 
describing the energy spectra of kaons and baryons ( antibar- 
yons )# All those are taken together but relative role of bar- 
yons and antibaryons increases with primary energy according 
to the expression 

F = 1 - .77(E o /205 )“*° 43 . ( 6 ) 

It is assumed that the average energy of kaons and baryons 
( antibaryons ) increases with energy faster than that of pio- 
ns and the relative inelasticity coefficient for production of 
those also increases with energy. The expression giving the 
ratio of the particles average energy and their inelasticity 
coefficient to those of pions are respectively 

b = ,45 + .253 InfF , R * 1.81 - 1 .5(E o /205 )~ #042 . (7) 

The Interactions of mesons with protons are also assumed 
to be described by the same expressions except that the total 
inelasticity coefficient is taken to be 0.6 instead of 0.5 
with total multiplicity the same and average energy increaesd 
by factor of 1.2 . 

Taking into account all these assumptions formula (l) is 
transformed into the following expressions for kaon and baryon 
plus antibaryon production in proton interaction and for pion 
production in pion interactions 

P(E)=:(E o /205) << q((E/b/E )(E /205)°S — , ( 8 ) 

E 

o/ k(E ,205) 

P(E)»(E o /205) q((-/l.2/E o )(E /205) ) 2 . ( 9 ) 

L E [_ 

Combining both formulae we can obtain expression for the spec- 
trum of kaons and baryons plus antibaryons produced in the me- 
son interactions. 

2 # APROXIMATE FORMULAE FOR INTERACTIONS WITH AIR NUCLEI 

For analysis of the cosmic ray propagation in the atmosp- 
here in actual fact, instead of the formulae for interactions 

of protons and mesons with protons 9 we have to use formulae a— 
ppropriate for interactions with air nuclei* Using the method 
outlined among others by Elias et al ( 1980 ) we have introduced 
simple corrections to the above given expressions to account 
for the fact that target is an air nucleus. 

Similarly as Elias et al ( 1980 ) we have assumed that an 
interaction of hadron with nucleus is an superposition of sub- 
sequent interactions of leading nucleon. The nucleon interact- 
ing with nucleus is assumed to undergo with certain probabili- 
ties one , two or more interactions with nucleons in the nucl- 
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eu$. The probabilities for the number of interactions used by 
us are given in the table X. The energy spectrum of the secon- 
daries and the inelasticity coefficient distribution the case 
of one interaction are the same as for interactions with prot- 
ons. When the number of interactions is higher those distribu- 
tions are obtained by summing the contribution from subsequent 
int erac t ions . 



, ( 10 ) 


where K. is the inelasticity coefficient in the i - th subseq- 
uent interaction. 

In similar way the inelasticity coefficient for p-nucleus 
interaction can be obtained. Overall spectra of secondary par- 
ticles can be easily obtained adding the contributions from 
various number of interactions with probabilities given in ta- 
ble I, The distribution of the inelasticity coefficient for pA 
collision can be obtained by Monte Carlo allowing the subsequ- 
ent inelasticities in p-nucleon interactions to fluctuate acc- 
ording to the assumed fluctuations of the inelasticity coeffi- 
cient in pp interactions. 


Table I 


Probability of an incident proton having n collision 
within a carbon nucleus 


n 

i 

2 

3 

4 

5 

P 

n 

.51 1 

.301 

.131 

.044 

• 

o 

U3 


Detailed calculations showed that the obtained energy spectra 
of the secondaries can be well represented by very simple mod- 
ification of the expression (l) and also the expressions ( 8 ) 
and ( 9 ). The modified formula ( 1 ) becomes 


The inelasticity coefficient distribution calculated by the 
outlined above method is consistent with an expression which 
we have developed ( Kempa and Wdowczyk 1985 ) on the basis of 
experimental data compiled by Jones (1983)* 


P(E W )«(FC) # 5 (E o / 205) q((E it /E o )(E o /205) FC°) 

FC = 1.32 


k(E .205) 


E 


. ( 11 ) 


rr 
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The expression is 

P(K) = 1.9exp(-2(l-K)), for K^.l and P(k)*VF , for K<,1 . 


Kempa and Wdowczyk ( 1985) gives also additional expression for 
the quasielastic interactions but those practically can be ne- 
glected here. Xn spite of the fact that they give almost 20/6 
of the total cross - section their contribution to the flux of 
secondaries is negligible. Their role practically amounts to 
an increase of the effective mean free path by a factor of 1.2 

Predictions based on the formula (ll) are compared with 
data of Elias et al (1980) in figure 1 where the distributions 
of pseudorapidities are given for various assumptions about 
relation of between p, . and p, . It seems that the best val- 
ue is obtained if we talce Tpp 

^ p tpA^ = 1,15 ^ p tpp^ * 

The above described picture of hadron-hadron interaction 
have been used in our calculations of cosmic ray propagation 
in the atmosphere presented elsewhere in this proceedings. 
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